Effect of hyperon-hyperon interaction on bulk viscosity and 
r-mode instability in neutron stars 

Debarati Chatterjee and Debades Bandyopadhyay 
Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata-700064, India 

Abstract 

We investigate the effect of hyperon matter including hyperon-hyperon interaction on bulk vis- 
cosity. Equations of state are constructed within the framework of a relativistic field theoretical 
model where baryon-baryon interaction is mediated by the exchange of scalar and vector mesons. 
Hyperon-hyperon interaction is also taken into account by the exchange of two strange mesons. 
This interaction results in a smaller maximum mass neutron star compared with the case without 
the interaction. The coefficient of bulk viscosity due to the non-leptonic weak process is determined 
by these equations of state. The interacting hyperon matter reduces the bulk viscosity coefficient 
in a neutron star interior compared with the no interaction case. The r-mode instability is more 
eff'ectively suppressed in hyperon-hyperon interaction case than that without the interaction. 

PACS numbers: 
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I. INTRODUCTION 



Neutron stars may undergo various instabilities which are associated with unstable modes 
of oscillation To understand oscillatory modes of a neutron star one has to study its 
stability properties. Each physical restoring force acting on a fluid element leads to a new 
family of oscillatory modes. Pulsation modes of neutron stars are classified according to 
different restoring forces. Here we are interested in Coriolis restored inertial r-modes. It was 
shown that r-modes were subiect to Chandrasekhar-Friedman-Schutz gravitational radiation 
instability in rapidly rotating neutron stars . This instability may play an 

important role in regulating spins of young neutron stars as well as old, accreting neutron 
stars in low mass x-ray binaries (LMXBs) and provides a plausible explanation for the 
absence of very fast rotating neutron stars in nature. This finds support from the study 
of eleven nuclear-powered millisecond pulsars whose spin frequencies are known^from burst 

Recently 



The 



oscillations has shown that the spin distribution has an upper limit at 730 Hz 0, |^ . 
observed fastest rotating neutron star is a radio pulsar with a spin frequency, 716 Hz 
r-mode instability could be a possible source for gravitational radiation |lO|. And, the 
detection of gravitational emission due to r-mode instability may shed light on the interior 
composition of a neutron star. 

It was argued that the r-mode instability could be effectively suppressed by bulk viscosity 
due to exotic matter in neutron star interior. Neutron star matter encompasses a wide 
range of densities, from the density of iron nucleus at the surface of the star to several 
times normal nuclear matter density in the core. Since the chemical potentials of nucleons 
and leptons increase rapidly with density in the interior of neutron stars, different kinds 
of exotic matter with large strangeness fraction such as, hyperon matter, Bose-Einstein 
condensed matter of antikaons and quark matter may appear there. The coefficient of bulk 
viscosity due to non-leptonic weak processes involving hyperons was calculated by several 



authors 



liilQ,Q,Q,Ii3 



It was noted that when the gravitational-radiation growth time 



scale was longer than the damping time scale due to hyperon bulk viscosity, the r-mode 
instability would be completely suppressed. Actually, a window of instability was found in 
the calculations of Ref . [3, 3| • Also, the impact of bulk viscosity due to unpaired and paired 
quark matter on the r-mode instability had been investigated extensively 

Earlier studies showed that hyperons might appear in neutron star matter around 2-3 



ipaired ana pair 
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normal nuclear matter density 



It was also noted that A hyperons appeared first 

.1 



in the system followed by S hyperons [20]. On the other hand S hyperons were excluded 
from the system because of a repulsive E hyperon potential depth at normal nuclear matter 



density 



2l|. Also, hyperon-hyperon interaction becomes important because the matter is 

jy the exchange of two 
23l |. It was found that 



hyperon-rich at high density regime. This interaction was accounted 
strange mesons, /o (975) (denoted hereafter as a*) and </> (1020) {22, 
the competition between attractive a* and repulsive fields made the overall equation of 
state stiffer at higher densities |j23|. In previous investigations of hyperon bulk viscosity, S 
hyperons appeared before A hyperons because an attractive E hyperon potential at normal 
nuclear matter density was considered. Those calculations also did not include hyperon- 
hyperon interaction. 

In this paper, we investigate the effect of hyperon-hyperon interaction on bulk viscosity 
coefficient, the corresponding damping timescale and the r-mode stability. Further, we use 
recent experimental data of hypernuclei to determine hyperon-scalar meson couplings in 
this calculation. This paper is organized in the following way. In Sec. II, we describe the 
model to calculate equation of state (EoS), bulk viscosity coefficient and the corresponding 
timescale. Parameters and results of our calculation are discussed in Sec. III. We summarise 
and conclude in Sec. IV. 



II. FORMALISM 

We adopt a relativistic field theoretical model to describe the /3-equilibrated and charge 
neutral hadronic matter. The constituents of matter are r?,,^, A, S+, S^, S'^, S^, H° of the 
baryon octet and electrons and muons. In this model, baryon-baryon interaction is mediated 
by the exchange of scalar and vector mesons and for hyperon-hyperon interaction, two addi- 
tional mesons - scalar a* and vector 0, |2^ are incorporated. Therefore the Lagrangian 
density for the hadronic phase is given by 

= "^^B ihfid^' -mB + gaBO- - gcoB^f^^^'' " gpBlp,^B ' f^) 
B 

+ ]^[d,ad^a-mla^)-U{a) 

+ \mlu^u^ - ip^, • p^^ + im>^ ■ + Cyy . (1) 
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The isospin multiplets for baryons B = N, A, S and S are represented by the Dirac spinor 
"^B with vacuum baryon mass ttlb, and isospin operator t^ and uj^^ and p^y are field strength 



tensors. The scalar self-interaction term 



2J is 



U{a) = ^g,a' + \g,a' . (2) 
The Lagrangian density for hyperon-hyperon interaction (£yy) is given by 

B 

+ \ {d,a*d^a* - ml,a*') - ^<P,^r' + ^mj^,^'^ . (3) 



We perform this calculation in the mean field approximation |25[. The mean values for 
corresponding meson fields are denoted by a, a*, uq, po3 and 0o- Therefore, we replace 
meson fields with their expectation values and meson field equations become 

dU 

+ ^ 

B 



mla = -— + ^gaBn% , (4) 



2 * 



^9a*Bn% , (5) 

B 

^gcvBriB , (6) 

B 

^1^0 = 9<j>BnB , (7) 



B 



\p02. = Y QpBhBnB ■ (8) 

B 

The scalar density and baryon number density are 

s _ 2 Js + 1 f^^B 2 r, 

UB = (2Jb + 1)^ , (10) 

where Fermi momentum is kpg, spin is Jb, and isospin projection is I^b- Effective mass and 
chemical potential of baryon B are m*^ = tub — gaB(^ — da'B^* and Hb = (kpg + m*^)^^'^ + 
9ujbuJo + g<pB<Po + hB9pBpo3, respectively. 



As the hadronic phase is charge neutral, the total charge density is 

Q = ^ Qbtib - He - = , (11) 

B 

where is the number density of baryon B, qs is the electric charge and Ug and are 
charge densities of electrons and muons respectively. 

In compact star interior, hypcrons maintain chemical equilibrium through weak processes. 
The generalised /9-decay processes may be written in the form Bi — > B2 + I + i^i and 
B2 + I — > Bi + i>i where Bi and B2 are baryons and 1 is a lepton. Therefore the generic 
equation for chemical equilibrium condition is 

IJ'i = hlJ-a - QilJ'e , (12) 

where yUg and //j are respectively the chemical potentials of neutrons, electrons and i-th 
baryon and bi and qi are baryon and electric charge of i-th baryon respectively. The above 
equation implies that there are two independent chemical potentials /i„ and corresponding 
to two conserved charges i.e. baryon number and electric charge. 
The total energy density is given by 

4/ {k^ + mfy/^e dk . (13) 
Jo 



And the pressure is 



^22^ 3 1 4 
-2^<^^ -3^2(7 --gso- 

1 2 *2|-'- 22,-'- 2 i2 , ^ 22 



1 2 Js + 1 dk 



+ qZ^^2 / Ci.2 I ™2Nl/2 • ^ 



(A:2 + mf)V2 
dk 

3 ^ 7r2 Jo (A;2 + m2)V2 
Bulk viscosity arises because pressure and density variations associated with the r-mode 
oscillation drive the system out of chemical equilibrium. Various microscopic reaction pro- 
cesses bring the system back to equilibrium. Weak interaction processes are most important 



in this case. As we are concerned about bulk viscosity coefficient in you ng n eutron stars 
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13| that 



where temperature is ~ 10^ - 10^'' K, it was shown by various authors 
non-leptonic processes involving hyperons might lead to a high value for the bulk viscosity 
coefficient. Though there may be several processes contributing to the bulk viscosity, only 
those reactions which might provide an upper limit on the bulk viscosity would be of interest. 
The relevant non-leptonic processes are 

n + p ^ p + A , (15) 
n + n ^ p + , (16) 
n + n ^ n + A . (17) 

All these reactions involve variation of neutron number density (n„) due to density pertur- 
bation on the system. Therefore, we consider neutron fraction as a primary variable. 

The real part of bulk viscosity coefficient is calculated in terms of relaxation times of 
microscopic processes 13, 13] 

1 + [iOT) 

where P is the pressure, r is the net microscopic relaxation time and and 70 are 'infinite' 
and 'zero' frequency adiabatic indices respectively. The factor 

dP dXn /I r>\ 

^--^' = -pd^.d^r ^''^ 
can be determined from the EoS. Here x„ = ^ gives the neutron fraction in the equilibrium 
state and Uh = ^^nB is the total baryon density. In the co- rotating frame, the angular 
velocity (a;) of (l,m) r-mode is related to angular velocity (fi) of a rotating neutron star as 



Here we consider the process in Eq. (15) because other processes may not contribute 
significantly to the calculation of bulk viscosity coefficient and are discussed later in details. 



The relaxation time (r) for the process is given by 
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2 

where pa is the Fermi momentum for A hyperons and < |Ma| > is the angle averaged 



matrix element squared given by Ref . jlO|, |l3| . Also, we have 

ni,6x„ 



ir vvuvf 

I3Q- 

Oinn - ttAn " ^nA + "AA , (21) 



where 5xn = Xn—Xn is the departure of neutron fraction from its thermodynamic equihbrium 
value Xn and = 1^ which is determined from the EoS. As soon as we know the 
relaxation time, we can calculate the bulk viscosity coefficient. 



The bulk viscosity d amp ing timescale (tb) due to the non-leptonic process involving 
hyperons is given by Q, Q 

1 _ 1 dE 

^ ~ ~2E1E ' ^^^^ 
where E is the energy of the perturbation as measured in the co-rotating frame of the fluid 

and is expressed as 

1 

E = -a^Q^R-^ / e{rydr . (23) 
2 Jo 

Here, a is the dimensionless amplitude of the r-mode, R is the radius of the star and e(r) 
is the energy density profile. The derivative of the co-rotating frame energy with respect to 
time is 

^ = -Att / C{r) <\V-6v\^> r^dr, (24) 
dt Jo 

where the angle average of the square of the hydrodynamic expansion j27[ is < | V ■ >= 
^^^(-^)^(l + 0.86(-^)^)(;^)^ and e is the mean energy density of a non-rotating star. Equa- 
tion (24) is the bulk viscosity contribution to the imaginary part of the frequency of the 
r-mode. Taking into account time scales for gravitational radiation (tgr) and modified Urea 
process {tjj) involving only nucleons, we can define the overall r-mode time scale (x^) as 

1 111 

- = + — + — . (25) 

Tr TgR Tb Tu 

Here is the imaginary part of the frequency of the r-mode. The gravitational radiation 



timescale is given by 



1310727r '•^ 



Q'' / e(r)rVr . (26) 
TGR 164025 Jo 

The time scale due to modified Urea process (t[/) involving only nucleons is calculated 
similarly as the hyperon bulk viscosity damping time scale but usin g th e following expression 



for bulk viscosity coefficient for modified Urea process given by j4, 28 1 



Cu = Qx lO-'h^T'^u^ . (27) 

Damping timescales are functions of angular velocity, neutron star mass and temperature. 
Therefore, solving = 0, we calculate the critical angular velocity above which the r-mode 
is unstable whereas it is stable below the critical angular velocity. 



III. RESULTS AND DISCUSSION 



In this calculation, we require the knowledge of coupling constants for baryons with cr, a;, 
p, cr* and mesons. Nucleon-meson coupling constants are obtained by reproducing nuclear 
matter saturation properties and taken from Ref. Here we consider the parameter set 
corresponding to the value of incompressibility of nuclear matter K = 240 MeV. Nucleon- 
meson coupling constants are listed in Table I. 

Hyperon- vector meson coupling constants are determined from SU(6) symmetry of the 



quark model 



The scalar a meson coupling to hyperons is calculated from the 



potential depth of a hyperon (Y) 



(28) 



in normal nuclear matter. The potential depth of A hyperons in normal nuclear matter 
[/^(tt-o) = -30 MeV is obtained from the analysis of ener gy l evels of A hypernuclei 



aer gy 



Recent S-hypernuclei data from various experiments [33|, |3^ give a relativistic potential of 
U'^ {no) = -18 MeV. However, the analysis of E~ atomic data implies a strong isoscalar 
repulsion for S~ hyperons in nuclear matter {2^. Also, recent S hypernuclei data indicate 
a repulsive S-nucleus potential depth |35|. Therefore, we adopt a repulsive potential depth 
of 30 MeV for S hyperons Q. 

The hyperon-cr* coupling constants are estimated by fitting them to a potential depth, 
^ (no) , for a hyperon (Y) in a hyperon (Y') matter at normal nuclear matter density 



obtained from double A hypernuclei data 



23|. This is given by 



up (no) 



2[/f)(no) = 2[/f)K) 



-40 . 



(29) 



Now we perform our calculation employing nucleon-meson and hyperon-meson coupling 
constants corresponding to K = 240 MeV case. We find that before the appearance of 
hyperons, the /3-equilibrated matter is composed of neutrons, protons, electrons and muons 
and their abundances increase with baryon density. Here, charge neutrality is maintained 
by protons, electrons and muons. Particle populations as a function of baryon density and 
without hyperon-hyperon interaction are displayed in Fig. 1. At 2.6no, A hyperons appear 
first in the system. It is observed that the density of A hyperon increases whereas the neutron 
density decreases with increasing baryon density for the growth of A hyperons occurs at the 
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expense of neutrons. Next S~ hyperons are populated around S.Ono. However, S hyperons 
do not appear in the system even at higher densities because of repulsive S-nuclear matter 
interaction. Further we note that threshold densities for the appearance of different hyperon 
species except H° hyperons are not changed appreciably due to hyperon-hyperon interaction 
and this is demonstrated in Fig. 2. 

Equations of state (pressure versus energy density) for different compositions are shown 
in Fig. 3. The bold solid curve implies nucleons-only matter whereas the dashed and solid 
curves correspond to hyperon matter with and without interaction, respectively. The ap- 
pearance of hyperons makes equations of state softer. On the other hand, the EoS with 
hyperon-hyperon interaction is initially softer compared with the case without the interac- 
tion, but it becomes stiffer later. This may be understood in the following way. Hyperon- 
hyperon interaction is mediated by a* and mesons. The additional attraction of a* field 
makes the EoS softer. The repulsive contribution of 4> field becomes dominant at higher 
densities. There is a competition between these two effects. Consequently, the EoS changes 
over from softer to stiffer at higher energy densities. 

Now we calculate the bulk viscosity coefficient using the EoS including hyperon-hyperon 
interaction. Already, we have noted that A hyperons are populated first in the system and 
heavier hyperons appear at higher densities. Therefore, we are interested in non-leptonic 
processes, Eqs. (15) and (17), involving A hyperons. However, it was noted that transition 
rate for the non-leptonic process in Eq. (17) was one order of magnitude greater than that of 
the process in Eq. (15) and led to a lower relaxation time in the former case Q]. Therefore, 
we only consider the process in Eq. (15). In this calculation, partial derivatives of pressure 
and chemical potentials with respect to particle number density are determined with the help 
of the EoS. In Fig. 4, the difference of fast and slow adiabatic indices (700 — 70) is plotted 
with normalised baryon density. The dashed curve represents the case with hyperon-hyperon 
interaction whereas the solid curve indicates the case without the interaction. The dashed 
curve rises to a maximum value and then drops below the solid curve. This can be attributed 
to the interplay of a* and fields and the crossover from a softer to a stiffer EoS. 

Next we compute the relaxation time as given by Eq. (20) for the process in Eq. (15). 
Effective masses, Fermi momenta, chemical potentials of baryons and partial derivatives en- 
tering into the expression of the relaxation time are obtained from the EoS. In the calculation 
of matrix element, we use the values of axial- vector coupling constants Qnp — ~ 

1.27 and g^x = 
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-0.72 obtained from /5-decay of baryons at rest. We exhibit the relaxation time as a function 
of normahsed baryon density in Fig. 5. The temperature is set to 10^° K in this case. The 
dashed and sohd curves denote cases with and without hyperon-hyperon interaction, respec- 
tively. We observe that the effect of the interaction on the relaxation time is not significant. 
Further, the behaviour of relaxation time with temperature for hyperon-hyperon interaction 
case is shown in Fig. 6. The relaxation time increases as the temperature is decreased. The 
effect of temperature on the relaxation time is substantial. 

The hyperon bulk viscosity coefficient is plotted with normalised baryon density for a 
temperature 10^° K in Fig. 7. For this calculation we obtain angular velocities (Q) at 
different central densities corresponding to a sequence of neutron stars using the rotating 



star model 



37| and equations of state with and without hyperon-hyperon interaction. 



We note that the factor ut is much much less than 1 over the range of baryon densities 
that we have shown here and is neglected in the bulk viscosity coefficient given by Eq. 
(18). The bulk viscosity coefficient in hyperon-hyperon interaction case represented by the 
dashed curve becomes smaller than that without the interaction denoted by the solid curve 
in the density regime which might occur in maximum mass neutron stars. The effect of the 
interaction enters into the bulk viscosity coefficient through the pressure term and adiabatic 
indices. The bulk viscosity coefficient increases appreciably with decreasing temperature as 
it is evident from Fig. 8. The large value of ( at 10^ K indicates that the viscous damping 
might be important in suppressing r-mode instability in neutron stars. It is to be noted here 
that we have neglected A-hyperon superfluidity in this calculation. Earlier calculations of 
bulk viscosity with hyperons involved hyperon superfluidity [lol ll^l . However, it has been 
reported that superfluid gaps for hyperons would be very small jsJSl- 

Now we discuss the effect of the EoS including hyperon-hyperon interaction on the sta- 
bility of r-mode. This mode becomes stable when the gravitational radiation timescale 
is greater than the damping timescales due to hyperon bulk viscosity as well as modified 
Urea bulk viscosity. The calculation of bulk viscosity timescales as given by Eqs. (22)-(24) 
needs the structure of a neutron star and its energy density profile. We have calculated 
the sequences of non-rotating neutron stars using Tolman-Oppenheimer-Volkoff equations 
with and without hyperon-hyperon interaction. The maximum neutron star masses with 
and without hyperon-hyperon interaction are 1.64 Mgoiar and 1.69 Mgoiar corresponding to 
central baryon densities 7.1 and 6.9 uq, respectively. Similarly, we have computed sequences 
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of rotating stars with and without hyperon-hyperon interaction using rotating neutron star 



model of Stergioulas 



36 



37|. In this case, maximum masses of rotating neutron stars 



with and without hyperon-hyperon interaction are 1.95 Mgoiar and 2.00 Mgoiar corresponding 
to central baryon densities 5.4 and 5.7 uq respectively. 

Next we perform the calculation of damping timescale due to hyperon bulk viscosity 
with and without hyperon-hyperon interaction for rotating neutron stars of mass 1.6 Mgoiar- 
This requires the knowledge of energy density profiles in rotating neutron stars. Earlier 
investigations adopted slowly rotating neutron stars and used the energy density profiles of 
non- rotating neutron stars 13]. Recently, Nayyar and Owen Q| incorporated the effects of 
rotation in the calculation of critical angular velocity using Hartle's slow rotation approx- 
imation. Here we also include the effects of rotation on damping time scale using energy 
density profiles of rotating neutron stars. A rotating neutron star has less central energy 
density than that of its non-rotating counterpart having the same gravitational mass because 
of the contribution from centrifugal force. Therefore, hyperon population will be sensitive to 
the rotation. Here we are dealing with a situation where the young neutron star is rotating 
very fast and radiating gravitational waves giving rise to the r-mode instability in it. As the 
neutron star emits gravitational radiation it spins down. Consequently, the central density 
in it increases. When the baryon density in the interior is such that hyperon thresholds are 
reached, hyperons will be produced abundantly. In this calculation, we consider 1.6 Mgoiar 
mass neutron stars with and without hyperon-hyperon interaction having central baryon 
densities 3.9 and 3.6 hq which are much above the threshold for A hyperons and rotating at 
angular velocities Qrot = 2952 and 3100 respectively. We calculate the hyperon bulk 
viscosity damping time scales using the energy density profiles of those stars. Their energy 
density profiles with and without hyperon-hyperon interaction are displayed in Fig. 9. The 
dashed curve including the interaction has larger energy densities than the case without the 
interaction. It is to be noted that the Keplerian angular velocities of neutron stars of mass 
1.6 Msoiar with and without the interaction are 5600 s~^ and corresponding central baryon 
densities are below the hyperon thresholds. The angular velocity for l=m=2 r-mode for a 
neutron star rotating with angular velocity is given by u; = In this calculation we 
obtain a set of values for u corresponding to Q ranging from to Qrot- Also, we express 
hyperon bulk viscosity C as a function of r using the knowledge of baryon density profiles. 
All these inputs are now used to calculate the damping timescale. Similarly, we calculate 
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the damping timescale due to modified Urea process for nucleons using tlie expression for 
the bulk viscosity coefficient Eq. (27). However, the bulk viscosity due to leptonic processes 
is several orders of magnitude lower than that of non-leptonic processes and has small effect 
on the damping of the mode. It is to be noted here that the effect of direct Urea process 
involving hyperons on bulk viscosity was investigated and it was found to have a small ef- 
fect on the damping of the r-mode As the gravitational radiation drives the r-mode 
unstable, it comes with a negative sign in Eq. (25). It is worth mentioning here that the 
energy density profiles of rotating stars used in this calculation do not differ much from their 
non-rotating counterparts. 

Now we determine critical angular velocities {^c) as a function of temperature for a 
neutron star of mass 1.6 Mgoiar by solving = as given by Eq. (25). The critical angular 
velocity is plotted with temperature in Fig. 10 for cases with (dashed curve) and without 
(solid curve) hyperon-hyperon interaction. We find that two curves do differ a little bit 
at and above 5 xlO^ K. This implies that gravitational radiation dominates in this region 
making the r-modes unstable. However, below this temperature damping timescales due to 
hyperon bulk viscosity starts dominating. The hyperon-hyperon interaction suppresses the 
instability more effectively below 5 x 10^ K. Consequently, the instability window shrinks. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the effect of exotic matter, in particular, hyperon matter including 
hyperon-hyperon interaction on bulk viscosity. Here we have constructed equations of state 
within the framework of a relativistic field theoretical model. As large number of hyperons 
may be produced in dense matter, hyperon-hyperon interaction becomes important and have 
been included in our calculation. This interaction is mediated by two strange mesons. Here, 
we use recent hypernuclei data which give rise to attractive potential depths for A and S 
hyperons and a repulsive potential depth for S hyperons in normal nuclear matter. Also, we 
exploit the knowledge of double A hypernuclei data to find the strength of hyperon-cr* meson 
coupling constant. Using these potential depths, we find that A hyperons appears first in the 
system followed by and S° hyperons. However, S hyperons do not appear because of the 
repulsive potential. Hyperon-hyperon interaction makes the EoS softer resulting in a smaller 
maximum mass neutron star than that without the interaction. Next we have computed 
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the bulk viscosity coefficient due to the non-leptonic weak process n + p ^ p + K and its 
influence on r-mode stability. We have used energy density profiles for rotating neutron 
stars to take into account the effects of rotation on hyperon populations. It is found that 
the gravitational radiation driven r-mode instability is more effectively suppressed due to the 
bulk viscosity coefficient in hyperon-hyperon interaction case compared with the situation 
without the interaction. 

Besides hyperon matter, there might be other forms of matter such as Bose-Einstein 
condensates of antikaons and quarks. It was shown in earlier calculations that Bose-Einstein 
condensation of antikaons might appear around 2-377-0 and delay the appearance of hyperons 
[4^. So far, there is no calculation of bulk viscosity in the antikaon condensed phase due 
to non-leptonic processes and how it competes with hyperon bulk viscosity coefficient. This 
problem is being investigated by us and will be reported in a future publication. 
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TABLE I: The nucleon-meson coupling constants in the GMl set are taken from Ref. [29]. The 
couphng constants arc obtained by reproducing the miclear matter properties of binding energy 
E/B = —16.3 MeV, baryon density no = 0.153 fm~^, asymmetry energy coefficient ciasy — 

32.5 

MeV along with incompressibility K = 240 MeV, and effective nucleon mass m*^/mN = 0.78. The 
hadronic masses are mjv = 938 MeV, = 550 MeV, m^, = 783 MeV, and nip = 770 MeV. All 
coupling constants are dimensionless, except g2 which is in fm~^. 
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Fig. 2. Same as in Fig. 1 but for hyperon-hyperon interaction case. 
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Fig. 3. The equation of state, pressure P vs energy density e , for nucleons-only matter (bold solid 
curve), hyperon matter with (dashed curve) and without (solid curve) hyperon-hyperon interaction 
are shown here. 
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Fig. 4. The difference of adiabatic indices, 700 — Jo, is shown as a function of normahsed baryon 
density with (dashed) and without (sohd curve) hyperon-hyperon interaction. 
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Fig. 5. Relaxation time is plotted with normalised baryon density for the non-leptonic process 
in Eq. (15) at a temperature 10^*^ K for hyperon matter with (dashed curve) and without (solid 
curve) hyperon-hyperon interaction. 



20 




Fig. 6. Same as in Fig. 5 but for different temperatures and with hyperon-hyperon interaction. 
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Fig. 7. Bulk viscosity coefficient is exhibited as a function of normalised baryon density for the 
process in Eq. (15) at a temperature 10^° K for hyperon matter with (dashed curve) and without 
(solid curve) hyperon-hyperon interaction. 
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Fig. 8. Same as in Fig. 5 but for different temperatures and hyperon-hyperon interaction. 



23 



1.25e+15 



K = 240MeV 




5 10 15 

r (km) 



Fig. 9. Energy density profile is shown with radial distance for rotating neutron stars of mass 
Mgoiar with (dashed curve) and without (solid curve) hyperon-hyperon interaction. 
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Fig. 10. Critical angular velocities for 1.6 Mgoiar neutron star are plotted as a function of 
temperature with (dashed curve) and without (solid curve) hyperon-hyperon interaction. 
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